Bicomplexes and Backlund transformations 

A. Dimakis 1 and F. Miiller-Hoissen 2 



1 Department of Financial and Management Engineering, University of the Aegean, 
31 Fostini Str., GR-82100 Chios, dimakis@aegean.gr 

2 Max-Planck-Institut fiir Stromungsforschung, Bunsenstrasse 10, D-37073 Gottingen, 

fmuelle@gwdg.de 

Abstract 

A bicomplex is a simple mathematical structure, in particular associated with com- 
pletely integrable models. The conditions defining a bicomplex are a special form of a 
parameter-dependent zero curvature condition. We generalize the concept of a Darboux 
matrix to bicomplexes and use it to derive Backlund transformations for several models. 
The method also works for Moyal-deformed equations with a corresponding deformed 
bicomplex. 

1 Introduction 

Backlund transformations arose in the 19th century in a differential geometric context [p]]. 
Despite of the fact that this is already a rather old subject, it is still the subject of quite a lot 
of recent articles. Indeed, in the case of many nonlinear equations, a Backlund transformation 
(BT) remains the only hope to construct sufficiently complicated exact solutions. 

The essence of the concept of a BT is basically the following (see M ^ ffl, EL for example). 
Suppose we have two (systems of) partial differential equations^ EQi[ui\ = and EQ 2 [u 2 ] = 0, 
depending on a variable U\ and its partial derivatives, respectively u 2 and its partial derivatives. 
A BT is then given by relations between the two variables and their partial derivatives which 
determine u 2 in terms of u\ such that EQ 2 [u 2 ] = if EQx[ui) = holds. So, given a solution 
of EQi = 0, it determines a corresponding solution of EQ 2 = 0. This will only be of help, of 
course, if the relations between u\ and u 2 are considerably simpler than the equations EQi = 
and EQ 2 = 0, and if some solutions of one of these equations are known. For instance, EQi = 
and EQ 2 = may be higher order partial differential equations and the BT only of first order. 
If EQi = EQ 2 , then such a transformation is called an auto- Backlund transformation. It can 
be used to generate new solutions of EQi = from given ones. How to find (useful) BTs? For 
most completely integrable models, ways to construct BTs are known. The existence of such a 
transformation is usually taken as a criterion for complete integrability. 

A large class of completely integrable equations in two (space-time) dimensions admits a 
zero curvature formulation 

U t -V x + [U,V} = (1.1) 



1 We may consider difference equations as well. 
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(see ||, in particular) with matrices U(u, A) and V(u, A) depending on a parameter A, besides 
the dependent variable u. Here, U t and V x denote the partial derivatives of U and V with 
respect to the coordinates t and x, respectively. (|1.1| ) is the compatibility condition of the 
linear system 

z x — U (u, A) z , Zt = V(u,X)z (1.2) 

with a (matrix) 'wave function' z. A special class of auto-BTs is then obtained via A-dependent 
transformations 

zf = Q(u, u', A) z (1.3) 

depending on two fields, u and u', which preserve the form of the linear system, so that 

z' x = U(u, A) z' , z' t = V(u', A) z . (1.4) 

Then, if u is a solution of the equation modelled by the zero curvature condition, u' is also a 
solution. Q has to satisfy the equations 

Q x = U(u')Q-QU(u), Q t = V{u')Q-QV{u) (1.5) 

and is called Darboux matrix J7|, It is not known, except for special examples, whether 
all auto-BTs of an integrable model, possessing a zero curvature formulation, can be recovered 
in this way. Many auto-BTs are known to be of this Darboux form, however (see [[F], [TT|, for 



example). In many cases, an ansatz for Q suffices which is linear in A ||. The important 
'dressing method' of Zakharov and Shabat [12j actually involves the construction of a Darboux 



matrix. 

In terms of the A-dependent 'covariant derivative' D = d — Udx — Vdt = d + A, where d 
is the exterior derivative on IR 2 (with coordinates x and t), ( |1.2|) and (|1.4j ) can be written as 
Dz = and D'z' = 0, respectively, and (|T7Tp becomes D 2 = which is F = dA + A A A = 0.0 
( |1.5|) can be rewritten as 

A'Q = QA-dQ (1.6) 

which, for an invertible Q, is the transformation law of a connection ('gauge potential') under 
a gauge transformation given by Q. This is equivalent to the covariance property 

D'Q = QD. (1.7) 

of the covariant derivative. This scheme can be generalized to hetero-BTs where Q relates two 
different zero curvature equations.^ 

In the special case where D depends linearly on A, it naturally splits into two linear operators 
which do not depend on A and we have an example of a bicomplex |14[ (and even a bi-differential 



2 This is sometimes called a Darboux transformation (see 0, for example), but should not be confused with 
the classical Darboux transformation | jlO| which, however, is indeed related to a Darboux matrix 

3 Without further conditions, F = is simply solved by A = g^ 1 dg with an invertible matrix g. If we require 
the dependence of A on A to be polynomial or rational, this results in nontrivial equations, however. 

4 Various relations between zero curvature conditions and BTs have also been discussed in @, |L 0, Oj, for 
example. 
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calculus |T5|] ). A bicomplex is an N -graded linear space (over R) M = © S>0 M S together with 
two linear maps V, D : M s — > M s+1 satisfying 



V 2 = 0, D 2 = 0, DV + VD = . (U 



Introducing an auxiliary real parameter A, these equations can be written as a generalized 
parameter-dependent zero curvature condition,^ 

(T?-AD) 2 = 0. (1.9) 

Typically, the maps V and D depend on a (set of) variable(s) u such that the bicomplex 
conditions hold if and only if u is a solution of a system of (e.g., partial differential) equations. 
The corresponding linear system is 

(D-AD)z = (1.10) 

(for z G M°), or a slight modification of it. The simple linear A-dependence of this linear system 
has to be contrasted with the, in general, nonlinear A-dependence of U and V in ( |1.2| ). Thus, at 
first sight, the bicomplex formulation looks like a severely restricted zero curvature condition.^ 
However, rewriting the bicomplex equations, if possible, in the form ( |1 . 1| ) and ( |1.2| ), in general 
results in a nonlinear A-dependence of the corresponding U and V. But in general it will not 
be possible to rewrite a given zero curvature formulation, showing a nonlinear A-dependence, 
in bicomplex form[], although it is possible that a bicomplex formulation exists for this model. 
Disregarding the A-dependence, the structure of the bicomplex equations is somewhat less 
restrictive than (|1.1| ). In a series of papers [0, [T3|, |TB|, [T7|, [T8|, [05], [SO], several known integrable 
models have been cast into the bicomplex form, some new models have been constructed, and 
it has been demonstrated, in particular, how conservation laws can be derived from it in the 
case of evolution-type equations. A bridge from bi-Hamiltonian systems to bicomplexes has 



been established in It should be stressed, however, that a bicomplex structure is much 

more general and does not presuppose the existence of a symplectic or Hamiltonian structure. 

The idea of a Darboux matrix is immediately carried over to bicomplexes. Let £>» = 
(Mj,Dj,Dj), i = 1,2, be two bicomplexes depending on variables U\ and u 2 , respectively. 
We propose the following definition. 

Definition. A Darboux- Backlund transformation (DBT) for the two bicomplexes £>i and B 2 is 
given by a A-dependent linear operator Q(X) : M\ — > M 2 such that^ 

(V 2 - A D 2 ) Q(A) = Q(A) (V 1 - A D x ) (1.11) 

for all A. This is an auto-Darboux-Bdcklund transformation if the two bicomplexes are associated 
with the same equation.^ <0 

5 There are examples in section 4 where the bicomplex maps V and D involve complex quantities but are not 
linear over C. The bicomplex equations can then only be cast into this form for real A. 

6 Actually many integrable models can be derived by reduction of the self-dual Yang-Mills equation which 
possesses a zero curvature formulation linear in A. 

7 Clearly, this applies in particular to zero curvature equations with a non-rational dependence on A (see ||, 
for example). 

8 The requirement of 'form-invariance' of U and V in the Darboux matrix formalism is here replaced by the 
requirement that Q preserves the linear A-dependence of the operator V — AD. 

9 If Q does not depend on A and is invertible, an auto-DBT reduces to an equivalence transformation of the 
two bicomplexes. 
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Since the bicomplex maps depend on the fields u\ and u 2 (which are solutions of the respec- 
tive field equations), so does Q as a consequence of It implies 



Q(A) (Pi - A DO 2 = (V 2 - AD 2 ) 2 Q(A) . (1.12) 

Hence, if U\ is a solution of the equation associated with B\, then (V 2 — AD 2 ) 2 Q(A) = 0. If 
Q(X) is invertible, an assumption which we will make in all our examples, this implies that 
u 2 is a solution of the equation associated with B 2 . If Q(A) is not invertible, this cannot be 
concluded, in general. 

Given three bicomplexes which are connected by DBTs Q 2 \ '■ Mi — > M 2 and Q 32 : M 2 — > M 3 , 
the composition Q 32 Q21 is also a DBT: 

(X? 3 - A D 3 ) Q 32 (A) Q al (A) = Q 32 (A) Q 2 i(A) (Pi - A D x ) . (1.13) 

Let us now consider four bicomplexes B iy i = 0,1, 2, 3, with DBTs Q w , Q 20 , Q31, Q 32 . Suppose 
we can solve the respective DBT conditions such that solutions Ui of EQi = can be expressed 
via the DBT Qij in terms of solutions of EQj = 0. The condition 

Q3lQl = Q32Q2 (1.14) 

is then in many cases strong enough to guarantee that a solution uq of EQ = is taken via 
Q31Q10 and also via Q 32 Q20 to the same solution u 3 of EQ 3 = (see Fig. 1). The condition 
( |1.14| ) is our formulation of the 'permutability theorem' (see [|1], |^|, for example). In the case 
of auto-DBTs, it leads to nonlinear superposition rules for solutions of the respective equation. 
Typically, a BT (obtained from a DBT) depends on some arbitrary constants. The condition 
(ECU) then enforces relations between the corresponding constants on the left (uq i— > u\ t— > m 3 ) 
and the right way (uq 1— > u 2 i— > 113) in the right diagram of Fig. 1. The usual formulation of a 
permutability theorem is that if these relations between the constants hold, then the two ways 
yield the same solution u 3 . 



M 3 u 3 




Fig. 1: 

A diagram of DBTs and a corresponding induced mapping of solutions in case 
of commutativity of the first diagram. 

The linear dependence on the parameter A in the bicomplex zero curvature formulation 
greatly simplifies the derivation of DBTs. Specializing to various models, one easily recovers 
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well-known BTs. For this purpose we make an ansatz 



N 



Q(\)=J2* k Q ik) (1-15) 



k=0 



with some N E N and QW not depending on A. The DBT condition (|1.11| ) then splits into the 
system of equations 



P 2 Q (0) 


- g {0) Pi 


= 






V 2 


-Q {k) V 1 


= D 2 Q (fe - 1) -Q^'Dj 


(fc = l,...,iV) 


(1- 


D 2 gw 


- Q {N) D 1 


= 0. 







We speak of a primary DBT when AT = 1, of a secondary DBT when N = 2, and so forth. 
The composition of iV primary DBTs is obviously at most an iV-ary DBT. Generically, it will 
be indeed an iV-ary DBT.0 There may be higher DBTs which are not obtained in this way, 
however. 

If two bicomplexes admit an invertible DBT, there is an equivalent DBT-problem with Q(X) 
acting on a single bicomplex space M and = I, the identity operator.^ Expressing the 
dependence of T>i on a solution U{ explicitly as ^[wj]; the first of equations ( |1.16| ) then requires 
X?i[mi] = T> 2 [u 2 ] for all solutions U\ and u 2 related via Q(X). Let us consider the special case 
where the equation under consideration admits a bicomplex B\ such that D\ does not depend 
on u. In this case we write Xq = 5 and obtain D 2 = 5. The DBT-problem then takes the form 
(5 — A D 2 [«2]) Q(X) = Q(X) (5 — A Di[tti]). Moreover, if we look for aufo-DBTs, the bicomplexes 
£>i = (M, S, Di) and B 2 = (M, 5, D 2 ) have to be equivalent. This means that the respective sets 
of bicomplex equations, which depend on u, must both be satisfied if u solves the equation for 
which £>i and B 2 are bicomplexes. If this holds for £>i, then an obvious way to achieve that it 
also holds for B 2 is to choose T) 2 [u] = Di[u] = D[u]. Since 5 is common to both bicomplexes, we 
should actually hardly expect Di[w] and ~D 2 [u] to differ in a non-trivial way, although exceptions 
may exist .0 This motivates to restrict the invertible auto-DBT condition for an equation, which 
possesses a bicomplex B = (M, 5, D) where 5 does not depend on u, to the form 

(5 - \ D[u 2 ]) Q(X) = Q(X) (5 - A D[m]) (1.17) 

with Q(°) = I. Then we are dealing with a single bicomplex only. This restricted form of the 
auto-DBT condition underlies all examples in section 4. 

How severe is the above restriction on the bicomplex B\ (and thus B)7 Splitting a given 
bicomplex map D as T> = 5 + B with a suitable operator 5 which is independent of a solution 
u and satisfies 5 2 = 0, the generalized curvature T = [6, B] + B 2 vanishes (see also section 
3). T generalizes the classical differential geometric formula for the curvature of a connection 

10 An exception appears in the Liouville example treated in section 2, where the composition of primary DBTs 
of the form ( |2.12D is again a primary DBT. 

n If Q(A) is invertible, then determines an equivalence transformation of bicomplexes. Introducing 
£>2 = {Q (0) y 1 V 2 Q^ and D' 2 = D 2 Q (0) , we have Q(A) {V 1 - AD X ) = (D' 2 - AD 2 ) Q(X) with Q(A) = 

(Q(°))- 1 g(A) : M x ->Mi. 

12 If u i— > u is a symmetry transformation of the equation for u, then we may also choose D 2 [m] = Di[u]. Note 
that a symmetry of the equation is not a symmetry (equivalence transformation) of an associated bicomplex, 
in general. The consequences of this observation have still to be explored. 



5 



1-form B if 5 is given by an exterior derivative on some manifold. In that case, it is well-known 
that a gauge transformation exists which transforms B to B' = so that T>' = 5. In the much 
more general bicomplex framework we do not have a corresponding theorem at hand, though 
an analogous result should be expected for relevant classes of bicomplexes (see section 5 for 
an example). This would mean that, if a bicomplex exists, then also a bicomplex with a map 
T> which is independent of the solution of the underlying equation. At least, this motivates a 
corresponding ansatz. In practice, however, it is often difficult enough to find any bicomplex 
for some equation and it is then hardly possible to find a concrete transformation to such 
a special bicomplex. Moreover, such a transformation may change the concrete form of the 
equation (cf section 5) which then possibly makes it difficult to apply corresponding results 
(e.g., concerning DBTs) to the original problem.^ Furthermore, one should keep in mind that 
interesting examples may exist for which the above special bicomplex form cannot be reached. 
Of course, in such a case the DBT method can still be applied, although the calculations will 
be more involved, in general. 

In flTTL [18], [19]] we constructed bicomplexes for various Moyal-deformed classical integrable 



models. Here, the ordinary product of functions is replaced by an associative, noncommutative 
*-product P2| . The definition of a bicomplex, and in particular ( |1.8| ), as well as our definition 
of a bicomplex DBT still applies. Also in this case a DBT provides us with a helpful solution 
generating technique, as will be demonstrated with an example in section 4. 

Section 2 treats the example of the Liouville equation and its discretization. In section 3 we 



elaborate our definition of a bicomplex DBT for a 'dressed' form of the bicomplex maps [|1J, [15 



and the restricted case where all maps act on the same bicomplex space M. Section 4 then 
shows how to recover auto-BTs for several well-known integrable models from a bicomplex 
DBT. Section 5 deals with the Harry Dym equation (see [23|, in particular) and section 6 
collects some conclusions. 



2 Example: Liouville bicomplex 

In many examples, the bicomplex space can be chosen as M = M° ® A n where A n = ©™ =0 A r 
is the exterior algebra of an n-dimensional real vector space with a basis £ r , r = 1, . . . , n, of 
A 1 . It is then sufficient to define the bicomplex maps T> and D on M° since via 

n n 

d( k..,a ii ---c)= E (D0u...v)r---e- (2.1) 

?l,...,ir — 1 il,---,*r— 1 

(and correspondingly for T>) they extend as linear maps to the whole of M.fj In the case of 
A2, we denote the two basis elements of A 1 as r, £. They satisfy £ 2 = = r 2 and £r = — r£. 

1. Liouville equation. Let M = C°°(1R 2 , 1R 2 ) <g> A2. Let x, y be coordinates on IR 2 and z x , z y the 

13 In particular, the nonlinear Schrodinger equation (see section 4.2.7) is known to be "gauge equivalent" to 
the Heisenberg magnet model (see Q, for example). This can be understood as an equivalence of bicomplexes 

14 In some examples, the £ r can be realized as differentials of coordinates on a manifold. In this way contact is 
made with ordinary zero curvature formulations (linear in A) of continuous integrable models. The generalization 
from the algebra of differential forms to an abstract Grassmann algebra is important, however, in order to treat 
relevant examples within the bicomplex framework. 
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corresponding partial derivatives of z G M°. We define 

T>z = z x £ + (cr + — I)z r , Dz = k e 2 ^ o-Z £ + (z y + 4> y a^z) r 
with a constant k, the 2x2 unit matrix I, and 



0-3 



f 
-1 



f 




(7_ 




f 



Then (M, T>, D) is a bicomplex for the Liouville equation 

D 2<t> 



'XL) 



fte 



(2.2) 



(2.3) 



(2.4) 



Let us now consider two such bicomplexes corresponding to two different choices Kj, % = 1,2, 
for k and corresponding solutions 0«. Then the DBT condition (|1.11|) reads [D, Q] = X (D 2 Q — 
QDi). With the ansatz ( 1.15Q for Q, this becomes 

[V, Q®] = , [V, Q {k) ] = D 2 - Q ( ^ 1} Dx (A; = 1, . . . , N), D 2 Q (A ° = Q (JV) D x . (2.5) 

Taking for general 2x2 matrices, the entries of which are functions of x and y, the first 
equation leads to 



Q<°> = a(y) I + b(y) a + 



(2.6) 



with functions a, b which do not depend on x. The k — 1 DBT condition in particular requires 
a y = 0, so that a must be a constant. Furthermore, it leads to 



1 



Q W = fl -~[b y + b((j) 1 + 2 ) y ] a 3 + c(y) a + + ra^ 
with an arbitrary function c(y) and functions /, r subject to 



fx 



(2.7) 

(2.8) 
(2.9) 

(2.10) 

where a does not depend on y. If a 7^ 0, the two equations ( |2.9| ) now reproduce a well-known 
BT for the Liouville equation: 



r = a (0 2 - 00, , r x = a [« 2 e 2<fe - Kl e 2 ^ 1 ] . 
The next equation in ( |2.5| ) requires in particular]^] r y = r (0 X + <p 2 ) y and thus 

r = a e 01+02 



a e 



</>l+</>2 



1 



h + <p2) x = ~ {K2 e^ 1 - «i e^-^ - a, 
a 



(2-11) 



where a has been absorbed via a rescaling of a (cf 0, for example). This is precisely obtained 
as the primary DBT condition with 



Q = / + Aae^ 2 a_ 



;2.i2) 



15 For N = 1 this is a consequence of the last equation in (2.5). For N > 1, it follows from the k — 2 equation. 
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The remaining freedom in the general solution for Q (where N > 1) can only restrict this BT. 
In particular, it cannot lead to different BTs. The case a = and b 7^ also does not lead to 
additional BTs. 

If Ki = 1 and k 2 = 0, ( |2.11| ) is a hetero-BT connecting solutions of the Liouville equation 



Yxy — e 2<?i with solutions of 4> xy = (which is the wave equation in light cone coordinates). 

Let us now explore the permutability conditions for N — 1. Using cPl = 0, ( |1.14| ) with 
( [2.12| ) becomes 



^ = 0o + ln( a "t'" aae r V (2.13) 

This allows to compute a solution 3 in a purely algebraic way from solutions 0o, 0i, 02, if the 
pairs (0o, 0i) and (0o, 02) satisfy ( [2.11|) . 



The fact that N > 1 does not lead to other BTs is a special feature of the Liouville example. 
Let us consider the case «j = 1, i — 1,2, in more detail. In general, the higher DBTs should 
be expected to be compositions of primary DBTs (see also the KdV example in section 4.1). 
Indeed, in the following we show that the composition of two Liouville BTs is again of the form 
( [2. 11| ). Differentiating the Liouville equation (|2.4| ) with respect to x, we find <p xxy = 2(p x (j) xy 



and an integration with respect to y leads to 

<f>xx = f(x)+<Px 2 (2.14) 

with integration 'constant' f(x). Multiplying the first part of the BT ( p.ll|) by e^ 1- ^ 2 , we 
obtain — (e* 1- ^ 2 )^ = ae* = a <pi xy . Integration of the last equation yields 

<f> 2 = <f> 1 -ln{k-a ( j> lx ) (2.15) 

with integration 'constant' k(x). The latter is determined by a via 

(k/a) x + (1 - k 2 )/a 2 = f (2.16) 

which follows with the help of the second BT part in ( [2.11| ). Let us now consider two BTs with 

02 = 01 - hl(A:21 - "21 0Lr) , 03 = 02 - hlfe ~ "32 M ■ (2.17) 

Eliminating 02 from the second equation with the help of the first, using ( |2.14| ) we find 

03 = 0i - ln(/c 3 i - a 3 i ix ) (2.18) 

where 

&31 = a 32 k 21 + k 32 (a 2 i + a 2 i x ) , k 31 = k 32 k 21 + a 32 (k 21x - a 2 i/i) ( 2 -19) 



solve (|2.16|) with f\ = <pi xx — <j>\ x . Hence, composition of Liouville BTs preserves their form. 



Remark. The infinitesimal version of the first of equations (|2.11| ) is 50 y = 15 a e 2 ^ +s ^ = fiae 2 ^ 



where $ denotes a variation. Using the Liouville equation with k — 1, this can be integrated 
with respect to y, so that 50 = <f) x fia — ~8k where fik(x) is an 'integration constant'. This is 
also obtained as the variation of ( |2.15| ) about a = and k — 1. Together with the variation of 
the Liouville equation, d(f) xy = 2e 2 ^^<f), we obtain ik = —\~8a x and thus 50 = 5a <p x + \Sa x . 
As a consequence, [5i, 5 2 ] = 5a 3 4> x + \ fia 3x with 5a 3 = 5ai ia 2x — 5a 2 ia.\ x . 



S 



2. Discrete Liouville equation. Let M = M° ® A 2 where M° is the set of maps Z 2 — > M. 2 . In 
terms of the shift operators (S x z)(x,y) = z(x + l,y) and (S y z)(x,y) = z(x,y + 1), we define^] 

Vz = (S x z -z)£+ (cr+SyZ -z)r, Dz = k e s ^ a_S x z £ + (e^"^ as S y z - z) r (2.20) 

with a constant k. Then T> 2 = and D 2 = identically, and DD + DD = turns out to be 
equivalent to 

e -<f>{x+l,y) e -<f>{x,y+l) _ e -tj>{x+l,y+\) e -<)>(x,y) _ K ^ (2.21) 

Introducing coordinates u = (x + y)/A, v = (x — y)/ A with A 2 = k and 4>(x, y) = —ip(u — A, v), 
this reads 



e 



tp(u,v- A) e f(u,v+ A) _ e tp(u-A,v) e <f(u+A,v) _ ^2 ^ 22) 



which is Hirota's discretization of the Liouville equation [33], 2jj|. Let us now explore the 
corresponding DBTs with = I. Then we have to solve the equation 

[D,Q (1) ] =D 2 -Di (2.23) 

which restricts to the form = ro_ with a function r. Furthermore, we obtain the 
following set of equations, 

r (x,y) = e <t >1< y x >y)- ( t> 1 ( x >y +l ) — e ^(x,y)-(t>2(x,y+l) (2.24) 

r(x,y + 1) = ^^'V+^-^^'y) _ e <i>i( x >v+i)-<l>i{ x :V) (2.25) 
r(x + l,y)-r(x,y) = k 2 e M x +i,y)+H(*,y) _ Kl e *i(*+i.iO+*i(*.v) . (2.26) 

Using <9 +y = 0(x, y + 1) — 0(x, y), the first equation can be written as 

r — e -9+y<i>l _ e -d +y <j>2 _ fgd+yfa _ gd+yfa^ e ~d+y (<t>l +<t>2) ^ 27) 

with the help of which we can convert ( |2.25| ) into 

S y r = re a+ » Wl+fe) . (2.28) 
This equation can be 'integrated' and yields 

r = a {x) e 4 * 1 ^ 2 (2.29) 
with an arbitrary function a(x). Together with ( [2 .241 ), it leads to the first part of the BT, 

e -S y 4>i e -<t>2 _ e -4>i e -s v 4>2 = a _ (2.30) 

The other BT part follows from ( 2.26| ): 

(S x a) e s ^ 1+5 ^ 2 - a e 4 "^ 4 " 2 = k 2 e Sx4,2+4>2 - k x e Sx4>1+ ^ . (2.31) 

For AT = 1, no additional equations arise from the remaining DBT conditions. For AT > 1, we 
could at most get restrictions on the above BT. In particular, no new BTs can show up, as in the 
case of the continuous Liouville equation. Taking ( |2.29| ) into account, the N = \ expression for 
Q is the same as for the continuum model. Hence we obtain the same permutability condition. 
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Here and in the following we also use the shift operators acting on functions via (S x 4>)(x,y) = 4>(x + l,y). 
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3 Darboux-Backlund transformations of dressed bicom- 
plexes 

All the examples presented in the next section possess a somewhat more specialized form of the 
bicomplex equations than what is allowed by the general formalism of section 1. For this class 
one can make some general observations which help to reduce the amount of calculations needed 
to elaborate the DBTs in concrete examples. The corresponding formalism is developed in this 
section. Clearly, this is of a more technical nature. In principle, given a bicomplex formulation 
of some equation, the formulas of section 1 are sufficient to work out the corresponding DBTs. 
In a given example, however, it may turn out to be very difficult to do it in a straight way. 

It is often convenient [TJ], |T3| to start with a trivialQ bicomplex and to use what we call 



'dressings' to construct non-trivial bicomplexes. Normally, such a 'deformation' of a trivial 
bicomplex results in too many independent equations, but there are two particular ways of 
introducing dressings (see (A) and (B) below) which keep some of the bicomplex equations 
identically satisfied. 

Let (M, S, d) be a trivial bicomplex and L the space of linear operator valued formsQ acting 
on M, i.e., for z G M and T G L we have Tz G M. On operators we define 

dT=[d,T], 5T=[5,T] (3.1) 

where [ , ] is the graded commutator.! 19 ] Then (L, 5, d) is again a bicomplex and, moreover, a bi- 
differential calculus^ A dressing of the bicomplex (M, S, d) is a new bicomplex (M, V, D) 
where 

Bz = dz + Az, Vz = 5z + Bz (3.2) 

with 'connection' 1-forms A,BeL. The conditions for (M, T>,D) to be a bicomplex impose 
the following conditions on A and B, 

F = dA + A 2 = 0, F = 5B + B 2 = 0, SA + dB + AB + BA — . (3.3) 

Introducing a real parameter A and 

d A = 5-Ad, A(X)=B -XA, (3.4) 

the three conditions (|3.3|) can be compactly written as a A-dependent zero curvature condition, 

F(X) = dxA(X) + A(X) 2 = (3.5) 

for all A. 

Now we consider two dressings with Ai, Bi, i = 1,2, and look for a DBT of the two resulting 
bicomplexes The DBT condition ( |1 . 1 1|) then becomes 

d\Q + A 2 (X) Q — Q Ai(X) = . (3.6) 

Trivial' in the sense that the corresponding bicomplex conditions are identically satisfied. 
18 The elements of M s are called s-forms. 

19 <5 and d are odd. For an even operator T, [d, T] = dT — Td. For an odd T, [d, T]=dT + Td. 
20 Besides S 2 = d 2 = d<J+ld = 0, 5 and d obey the Leibniz rule, i.e., the graded product rule of differentiation. 
21 Note that in this framework the linear spaces Mi of the two bicomplexes are taken to be the same. If Q is 
invertible, this is not a restriction as pointed out in section 1. 
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In terms of A and B, this reads 



6Q + B 2 Q-QB 1 = XD 21 Q (3.7) 

where 

D 21 g = dQ + A 2 Q-QA 1 . (3.8) 

In the following we confine our considerations to the case where B% = = B 2 ^\ Then the 
conditions ( |3.3| ) for the two bicomplexes reduce to 

Fi = dAi + A\ = , 5A i = (3.9) 



8Q = \t) 21 Q. (3.10) 



and (|3.7p becomes 



This equation has to be solved in order to determine the DBTs of a dressed bicomplex. Using 
the ansatz (|1.15|) , the DBT condition splits into the following set of equations, 

6Q<® = 0, 6QW =D 21 Q( fc - 1 ) (k = l,...,N), D 21 QW = 0. (3.11) 

Now we assume that Q is invertible withQ = I and consider in more detail the case of 
a primary DBT with 

Q = I + XR (3.12) 
where R = does not depend on A. ( |3.11| ) then reduces to 



SR = D 21 I = A 2 - Ax , f) 21 R = . (3. 13) 

In fact, in all the examples which we explored so far, it turned out to be sufficient to consider 
such a primary invertible DBT in order to recover well-known BTs. We can somewhat simplify 
the last set of equations as follows, using two obvious ways to reduce the set of bicomplex 
equations 0.Q 

(A) Let Ai = g~ x dgi where gi : M° — > M° are invertible operators (e.g., matrices) not depend- 
ing on A. This solves the first of the bicomplex equations (|3~9"|). The second of equations ( |3. 13| ) 
is then equivalent to d(g 2 Rg^ 1 ) = which can be converted to da = by setting R = g^ 1 a g\ 
with a € The first of equations ( |3.13| ) now becomes 

$(92 la 9i) = 92^92 ~ 9i l dgi = g^dig^ 1 ) gi . (3.14) 



22 Actually, in general, this can be achieved by separate gauge transformations and coordinate transformations 
of the two bicomplexes, since T% = 0. See section 5 for an example. 



23 This choice for Q^> trivially solves the first of equations ( 3.11 ). In the case of auto-DBTs, arguments for 



this choice (under certain assumptions) have been given in section 1. 

2 The following is specific to bicomplexes and has no analogue in case of a zero curvature condition non-linear 
in A. Though such a condition can also be solved as in (A) by writing the A-dependent gauge potential as a 
"pure gauge" , the corresponding g then depends on A. 

25 For arbitrary N, this holds with R replaced by Q^ N K 
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(B) Let Ai = 5wi with Wi G L° not depending on A. This solves the second of the bicomplex 
equations (|3.9|). From the first equation of ( |3.13| ) we obtain 



R = W2 — w\ + T 



with 5T = 0. Then the second equation in (|3.13| ) becomes 

$2 - $1 + cLT + 5 (w 2 T - Tw\ + ^(w 2 - wi) 2 + ~[wi,w 2 



where 



$i = dwi + - [Swi,Wi] 



are 5-potentials of the curvatures F^, i.e., 

5& 



Fa . 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



Let us assume that the first 5-cohomology is trivial. The last equation together with Fi = 
then implies 



(3.19) 



Furthermore, 6dT = leads to dT = Sb with b G L°. Then ( |3.16|) can be integrated and leads 



\& 2 - ^1 + w 2 T - Twi + ~(w 2 - wi) 2 + ~[wx,w 2 ] +6 = 0. 



(3.20) 



In concrete examples, it is often simpler to work out directly the second of equations ( |3.13|) , 
however. 

For a DBT with Q of the form ( |1 . 1 5| ) where = I, the permutability condition ( |1.14j ) 
results in the following system of equations: 



Q31 + Q10 Q32 Q20 





k-i 



V31 V10 V32 V20 ~~ 2-~* v 2 



l V31 V10 V32 V20 



'20 

m=l 

fc = jV + l,...,2iV. 



m+n— fc 
Km,n<W 



For a primary DBT with Q = I + A R, this reduces to 

-R31 + -Rio = R32 + -R20 
R31 Rio — R32 R20 ■ 

In case (A), we have Rij = g^aij gj, so that ( |3.25|) becomes 

0-31 a 10 — &32 ^20 • 

In case (B), we have R^ = 0j — <pj + Tij. Then ( |3.24j ) reduces to 

T31 + T 10 = T32 + T 2 q • 



(3.21) 
N (3.22) 

(3.23) 

(3.24) 
(3.25) 

(3.26) 

(3.27) 



21. 



An 'integration constant' c with Sc — can be absorbed by a redefinition of b. 
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4 Bicomplexes and auto-Darboux-Backlund transforma- 
tions for various integrable models 

In the following, we elaborate auto-DBTs for various integrable models. All the examples of 
non-trivial bicomplexes in this section are of the form (M, 5, D) where (M, 5, d) is a trivial 
bicomplex and D has the decomposed 'dressed' form D = d + A as considered in the previous 
section. Assuming Q to be invertible, in section 1 we have motivated a restriction of the auto- 
DBT condition to the form ( 11.17 ) with = I. This is the basis for the following calculations. 



As in section 2, A2 denotes the exterior algebra of a two-dimensional real vector space. 

4.1 KdV and related equations 

1. KdV equation, primary DBT. Let M = C°°(R 2 ,R) <g> A 2 with 

6z = -3z xx t + z x £, dz = (z t + 4z xxx ) t - z xx £ (4.1) 
for z G M°. Dressed with the gauge potential 1-form 

A = 5w = -3(w xx + 2w x d x )T + w x £ (4.2) 
we get a bicomplex for the KdV equation 

u t + u xxx - 6uu x = (4.3) 



where u = w x (see also |16||). Here we choose w & L as a function (which acts by multiplica- 



tion). Looking for a primary DBT, we have ( |3.15|) with T x = 0. Furthermore 



I) 21 R = {R t + 4R xxx + 12R xx d x + 12R x dl-3[(w 2xx + 2w 2x d x )R 

-R (w lxx + 2w lx d x )}} t - (R xx + 2R X d x -w 2x R + Rw lx ) £ (4.4) 

has to vanish. Using ( |3.15| ), the £-part leads to 

(w 2 - wi) xx + 2(w 2 - wi) x d x - w 2x (w 2 -w 1 + T) + (w 2 - w 1 + T) w lx = . (4.5) 

In particular, this implies T = 2d x + (3 with a function (3{t) which, however, can be absorbed 
via a redefinition of w\ (which leaves the KdV equation invariant). Hence 

R = W2 - wx + 2d x (4.6) 

and ( |4.5| ) leads to the BT part 

(w 1 + w 2 ) x = 2a + ^(w 2 -w 1 ) 2 (4.7) 

with an integration 'constant' a which is an arbitrary function of t. The vanishing of the r-part 
of (0) together with (g77|) yields^ 

(w 2 - wi) t + (w 2 - Wi) xxx - 3(w 2 - wt) x (wi + w 2 ) x = (4.8) 



27 This is also obtained by integration of the difference of the two KdV equations for u\ = w\ x and ui = wi x 
with vanishing integration 'constant'. 
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which is the second BT part for the KdV equation (see H, p. 113, for example, and also 



26|, p7|). Introducing 

r = w 2 — W\ (4.9) 

the BT can also be written as 

w lx = a-^r x + ^r 2 , r t + [r xx - ^r 3 - 6ar] x = . (4.10) 
The last equation has the form of a conservation law and can be integrated once if we set 

r = -2(\n X ) x - (4.11) 

Then, in terms of \i the BT reads 

. Xxx Xt r Xxxx n Xxx X f a m\ 

w lx = a -\ , — = 6a ho h 7 — (4-12) 

A.X Ax A A.X 

with a 'constant of integration' "f(t), assuming \x 7^ 0. 

The following will be needed below for our discussion of the secondary DBT. We consider 
two BTs, 

(w 1 + w 2 ) x = 2oti + -(wx - w 2 ) 2 , (w 2 + w 3 ) x = 2a 2 + ^{w 2 - w 3 ) 2 . (4.13) 
The difference of both equations is 

(wt - w 3 ) x = 2(«i - a 2 ) + ^(wt - w 3 )(wx +w 3 - 2w 2 ) . (4.14) 
Solving for w 2 yields 

1 s x 2(«i — a 2 ) 

w 2 = -{w 1 + w 3 ) 1 , s = Wi-w 3 . (4.15) 

2 s s 

Inserting this expression into the sum of the two equations Q4.13j ), we obtain 

/ 1 \ s xx 1 fs x \ 2 1 2 2(a 1 -a 2 ) 2 

(wx + w 3 ) x = an + a 2 H (— J +-s H . (4.16) 

The complementary parts of the above two BTs are 



(wx - w 2 ) t + {wt - w 2 ) xxx - 3 {w{ x - wl x ) = 0, (4.17) 

(w 2 - w 3 ) t + (w 2 - w 3 ) xxx - 3 (wl x - wl x ) = . (4.18) 

Adding these two equations leads to 

st + s xxx - 3s x (wx + w 3 ) x = . (4.19) 

Using (|4.16| ) to eliminate (wx + w 3 ) x , we obtain 

1 1 3s 2 

s t + [s xx - 3(ax + a 2 )s- -s 3 + 6{ax + a 2 f- - -^] x = (4.20) 

14 



which, setting s = — 2(lnx) x . and integrating once with integration 'constant' 2e(i), becomes 

Xt of , \ c 3(«i - a 2 ) 2 / x V , X /, n 

3(0!! + or 2 ) - o x x ^ ( — ) +e— • (4.21) 

6r 



Xa; ^ V / X; 



Here, S x denotes the Schwarzian derivative 

XxXX 3 / XxX \ c~\c~\\ 

S ' X = — -2\TJ ■ (422) 

Correspondingly, ( [4.16 ) can be rewritten in the form 

«i + tt2 . Xxxx (Xxx\ 2 . (ui~a 2 ) 2 ( xY (A OQ x 

From these expressions one recovers for a\ = a 2 and e = a BT found by Galas J^J. The 
latter is therefore just the composition of two 'elementary' BTs. 

2. KdV equation, permutability. With (|4.6j ), the permutability condition (|3.24|) is identically 
satisfied and (|3.25| ) leads to 

2 (wi — W 2 ) x 

w 3 = -w + w 1 + w 2 - ^ . (4.24) 

Wi - w 2 

With the help of ( |4.7|) , for the pairs w , W\ with a\ and w , w 2 with a 2 , the last equation can 
be written as 

W3 = Wo - 4 ai ~ a2 . (4.25) 
Wi - w 2 



3. KdV equation, secondary DBT. We consider again the above bicomplex associated with the 
KdV equation. But now we turn to the secondary DBT. Again, we have = w 2 — Wi + T 
with T x = 0. For N = 2, ( |3 . 1 1[ ) then requires 



d(w 2 — w 1 + T) + (5w 2 ) w 2 + w\ 5wi + 5(-wiw 2 + w 2 T - Tw\) 
~5w 1 -{5w 2 )Q { - 2) . 



The £-part of ( |4.26|) can be integrated and leads to 



Q 



(2) 



(wi - w 2 ) x + -(wi - w 2 ) 2 + 2(wi - w 2 ) d x + w 2 T - Twi + p 



where p x = 0. Inserted in the £-part of (|4.27|) , which is 

Qi 2 j + q (2) Wlx - w 2x q (2) + 2 q(?) d x = o 



(4.26) 
(4.27) 



(4.28) 



(4.29) 



this enforces T = 4d x + Pit). The function (3{t) can be absorbed by a redefinition of w\. 
Furthermore, we obtain 



q(2) = _ rx + _ r 2 _ Awix + 2r d x + Adl + Aa 



(4.30) 
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with r = w 2 — W\ and an arbitrary function a(t), and finally 



(tin + w 2 ) x = 2a + ^ - \ (?*) 2 + \r 2 - ^ (4.31) 
r 2 V r / 8 



with an integration 'constant' j(t). The r-part of (|4.26| ) is now evaluated to 

r t + r xxx - 3r x (r x + 2w lx ) = . (4.32) 
Elimination of w\ x with the help of ( 4.31|) leads to 

r t + [r xx - - Qar - K 3 - 2A\ = . (4.33) 
2 r 8 r 

Setting r = — 2(\nx)x, this equation can be integrated and rewritten as 

^ = Qa-S xX + Ql(—) +e— (4.34) 

Xx \ Xx J Xx 

with an arbitrary function e(i) and the Schwarzian derivative defined in ( [4.22 ). In terms of \i 
( 4.31] ) takes the form 

/ \ 2 / \ 2 

Xxxa; / Xxa; \ IX 



. /\xxx i /\xx \ I A. I /i r\ 

Wi ' =a+ wr\^j -Atj ■ (435) 

Comparison with our previous results shows that this secondary DBT for the KdV bicomplex 
is just the composition of two primary DBTs. If 7 < 0, this is evident. If 7 > we set 
«i = a + 2^7 and a 2 = a — i^fy. Although in this case the elementary BTs (with a,\ and a 2 , 
respectively) do not produce real solutions from real solutions, in general, their composition 
does. 

4. Modified KdV equation. Supplied with the product (a, b)(a', b') = (aa f , bb'), the vector space 
1R 2 becomes a commutative ring with unit (1, 1) which we denote as 2 M [29|. It is a realization 
of the abstract commutative ring generated by a unit 1 and another element e satisfying e 2 = 1. 
Here we have e = (1, — 1). The relevance of 2 M for the mKdV equation stems from the following 
observation (see also PD[ ). Let u be a field with values in 2 M which satisfies the KdV equation 

u t + u xxx - Quu x = (4.36) 

and let 

u = v 2 + v x e (4.37) 

with a real valued field v. Now 

u t + u xxx - 6uu x = (2v + e d x )(v t + v xxx - 6v 2 v x ) (4.38) 
shows that the field v satisfies the mKdV equation 

v t + v xxx - 6v 2 v x = . (4.39) 
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U = [U,U 



In fact, ( |4.37| ) is the famous Miura transformation and its 'conjugate' since with 
(v 2 , v 2 ) + (v x , —v x ) we obtain 

u + = v x + v 2 , u~ = -v x + v 2 . (4.40) 

Consequently, the two KdV equations for u^ 1 are equivalent to the above mKdV equation. 
Hence, in order to find an auto-BT for the mKdV equation, we simply have to extend our KdV 
treatment to fields with values in 2 M, though we have to take care of the fact that 2 R is not a 
division ring (since (1, 0)(0, 1) = (1 + e)(l — e) = 0). Introducing w = v + ve with v x = v 2 , we 
have w x = u and we can directly generalize the KdV auto-BT: 

(wi + w 2 ) x -2a--{w 1 - w 2 ) 2 = (4.41) 
(wi - w 2 )t + (wi - 102)^ - 3{wi - w 2 ) x (wi + w 2 ) x = (4.42) 

(where (3{t) has been absorbed in w±). Note that the integration 'constant' a(t) is now an 
element of 2 R. We decompose (|4.41|) with a = —k 2 + be to obtain the two equations 

(v 1 + v 2 ) x = -2fc 2 + i[(t>! - z) 2 ) 2 + - i; 2 ) 2 ] (4.43) 
(v l + v 2 ) x = 2b + (vi-v 2 )(v 1 -v 2 ). (4.44) 

Using v x = v 2 in ( |4.43 ), we obtain 

( Vl + V2 f = -4k 2 + (ik - v 2 f . (4.45) 

Applying d x to this equation and comparing the result with ( |4.44| ), we find 6 = 0. Solving the 
last equation for v\ — v 2 and inserting this expression in (|4.44 ), leads to the first part of the 
auto-BT,Q 

(vi + v 2 ) x = ±( Vl - v 2 ) ^Ak 2 + ( Vl + v 2 ) 2 . (4.46) 

Decomposition of ( 4.42 ) leads to 

(vi - v 2 )t + (vi ~ v 2 )xxx - 3 (vi - v 2 ) x {vi + v 2 ) x - 3 (vx - v 2 ) x (vi + v 2 ) x = (4.47) 
{v\ - v 2 ) t + {vx - v 2 ) xxx - 3(«i - v 2 ) x (vi + v 2 ) x - 3 (vt - v 2 ) x {vi + v 2 ) x = 0. (4.48) 

Using v x = v 2 in the second equation produces the second part of the mKdV auto-BT, 

(vi - v 2 ) t + [(«! - v 2 ) xx -2v\ + 2 v% = (4.49) 

which is the difference of the two mKdV equations for v\ and v 2 . 

5. mKdV equation, permutability . In the framework of the previous subsection, the KdV 
permutability condition takes the form 

(w 2 - wi)(it>i + w 2 - w - w 3 ) = 2 (w 2 - wi) x (4.50) 



28 With v = q x , this reads (qi + q2) x = 2fc sinh(qi — 52) where an integration constant has been absorbed by 
a redefinition of the qi. 
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(cf (fO§) and, by use of fl£3ip , 

(w 3 - w ) (w 2 -w 1 )=A {k\ - k\) (4.51) 

from which we obtain, by decomposition, 

{v 3 -v )(v 2 -vi) + {V2~ v t )(v 3 -v ) = (4.52) 
(v 3 -v )(.V2-v 1 ) + (v 3 -vo)(v 2 -v 1 ) = A(k 2 2 -kl) (4.53) 

and thus the following superposition formula for mKdV solutions, 

4{k 2 k 2 ){v2 _ Vi) 

v 3 = v + 7 rr, 77 ^ . 4.54 

(v 2 -v x y - (v 2 -vi) 2 

6. ncKdV equation. We choose the bicomplex maps and the dressing as for the KdV equation, 
so that ( J4.1|) and Q4.2| ) hold. But now we take u to be a map from R 2 into some noncommu- 
tative associative algebra with product * for which d t and d x are derivations. Then we have a 
bicomplex iff dA + A* A = which is equivalent to the noncommutative KdV equation (ncKdV) 

u t + u xxx - 3(u * u x + u x * u) = (4.55) 



where u = w x \ V7\ . The corresponding potential ncKdV equation is 

Wt + w xxx — 3w x * w x = . (4.56) 



The conditions for a primary DBT with Q of the form (|3.12j ) are 



5R = A 2 - A x , dR + A 2 * R - R * A x = . (4.57) 

As in the ordinary KdV case, the first equation is solved by R = w 2 — ui\ + T with T x = 0. The 
second equation implies T = 2d x + (3{t). Again, the function (3{t) expresses the freedom in the 
choice of the potential for u and can be set to zero. Furthermore, we obtain 

(wi + w 2 ) xx -wi* w lx - w 2x *w 2 + (w 2 * wi) x = (4.58) 
(w 2 - w x ) t + (w 2 - w x ) xxx - 3(w 2x * w 2x - w lx * Wi x ) = 0. (4.59) 



In the commutative case, ( |4.58| ) can be integrated (which introduces a parameter in the BT). 



This is not so in the noncommutative case. But we still have a BT which is no longer symmetric 
in wi and w 2 , however. 

7. ncKdV equation, permutability. The permutability conditions are reduced to -R31 * R 10 = 
R32 * -R20 This yields 

w 3 * (w\ — w 2 ) + (w\ — w 2 ) * w = wi * ui\ — w 2 * w 2 — 2(wi — w 2 ) x (4.60) 

where each of the pairs (wo,wi), (wo,w 2 ), (wi,w 3 ), (w 2 ,w 3 ) has to satisfy the BT equations 
and ( f4.59p . Let us now specify the ^-product as the Moyal product 



/ *h = moe® p/2 {f ®h) (4.61) 
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for smooth functions /, h, where m(/ <g> h) = fh, P = d t ® d x — d x ® d t , and d is a deformation 
parameter. As an example, let w% = 0, so that 

(wi — W2) * wo = W\ * w\ — W2 * W2 — 2(wi — W2) x ■ (4.62) 

For Wi and W2 we choose the 1-soliton solutions 

wi = -2 tanh(a; - At) , w 2 = -A coth(2x - 32t) (4.63) 

(see also 0, p. 116). (wi, W3) and (wi, W3) indeed satisfy (|4.58|) and (|4.59|) . Then the *-products 
on the right side of ( |4.62| ) reduce to ordinary products and we obtain 

f(x,t)*w = g(x,t) (4.64) 



with 



This implies 



and 



g(x,t) = 4tanh 2 (x-4t) - 16 coth 2 (2x - 32t) -2f x (4.65) 
f(x,t) = -2 tanh(x - At) + 4 coth(2a; - 32t) . (4.66) 



= * wo) = f * ^ + ~(/ t * w 0x -f x * wot) (4.67) 



d 2 w dw 0x dwot, 1,, , , , , . 

■> * ~^2~ = ~U* * ~ * -Qfl-\ ~ 1\f*t * W 0xx ~ 2ftx * Woxt + Jxx * Wott) • (4.68) 

For vanishing deformation parameter the solution of the permutability conditions is the 
2-soliton solution 

Wo = -6/ [2 coth(2 x-32t)- tanh(x -At)] (4.69) 
with corresponding KdV solution 

3 + cosh(4x-64t)+4cosh(2s-8t) , 

uo = W 0x = -12 / 4.70 

(cosh(3 x - 36 1) + 3 cosh(a; - 28 1)) 2 

(cf [[|, p. 116). The noncommutative solution is then 

1 





Wo 


dw \ 


1 


d$ J 




d 2 w " 




d$ 2 , 


) -7 



-<? 2 

2 



(4.71) 



where 

^1 : = (^) =-sr(/t^ox-/ a! Wbt) (4-72) 
^ ==|^| = -7(/- ^ - /. W lt) - ±(/ a Wi. - 2/, ^ + /. ^) ,4.73) 



N 2 



and so forth. In the case under consideration, we obtain W\ = and 

- 331776 ( cosh ( 3 g~36t)-3 cosh(x - 28 1)) (sinh(3 x - 36 i) + sinh(a; - 28 1)) 2 
2x ~ (cosh(3x-36t) + 3 cosh(x - 28 1)) 5 V± ''' ±) 

which is precisely the expression for the second order ncKdV correction u 2 to the classical 
2-soliton solution ( |4.70[ ), obtained in [I7| in a different way. 



19 



4.2 Further examples 



1. Sine-Gordon equation. Let M = C°°(M 2 ,C) ® A 2 and 



= z x £ + z)t, dz=~^(z-z)£ + z y T 



(4.75) 



for z G M°, where a bar indicates complex conjugation. Dressing d with A = g _1 dg, where 
g = e - * 9 ^ 2 with a real function 0, we obtain the map Dz = e^/ 2 d(e~ l ^^ 2 z). Then (M, 5, D) is a 
bicomplex associated with the sine-Gordon equation <p xy = sin [|14[] . Following scheme (A) of 
section 3, in order to calculate the primary DBT, we have to evaluate ( 3.14j) with Oj = e~ %( ^^ 2 
and some operator a. The latter has to satisfy da = which means Hz = az and a y — 0. First 
we calculate the right side of ( |3.14| ): 



g^dfagl 1 ) 9i z = \{e^ - e^)zi - |(0 2 - <f> x ) v zr 



(4.76) 



In order for this to be consistent (for all z) with the left side of ( |3.14| ), the operator a must 
have the form az = az with a real function a(x). Then 



~5{g 2 l a gi ) z = [y (0i + (j> 2 ) x + a x ] e^+^z^- la sin 



+ 



Z T 



(4.77) 



and p. 14 ) results in the following two equations, 

2 . 02-0i , _ . a x 
— + 2i — , 



(0i + <\>2)x = - sin - — — 
a 2 



a 



2a sin 



+ 02 



(4.78) 



Since we consider only real sine-Gordon solutions, we have to set a x = 0. Hence a has to be a 
real constant. Now we recover a famous auto-BT of the sine-Gordon equation ^ 0, [3], 31 . 
Since 



(4.79) 



the permutability condition (|3.26|) is satisfied with a±o = a^ 2 = on and a 2 o = a%\ = a 2 . The 

(4.80) 



remaining permutability condition ( |3.24| ) requires 

a x e*^"^)/ 2 + a 2 e ^+<«/ 2 _ a2 e i(*>+fc)/2 _ ai jfo+te)/* = Q 

which is equivalent to Bianchi's 'permutability theorem' for the sine-Gordon equation: 

«i + a 2 



+ 4 arctan 



tan 



a x — a 2 



4 



(4.81) 



This determines a solution 03, if 0i and 02 are obtained from 0o via the BT, i.e. the pairs 
(0o, 0i) and (0o, 2 ) have to satisfy (|4.78|) .H 



29 The reader should be aware of a problem of notation in this section. A BT like ( 4.78 ) is written in terms 
of solutions <pi and <p2- But these are not the solutions (f>i and </>2 which appear in the permutability relations. 
See also section 1. 
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2. An equation related to the sine-Gordon equation. Let us again consider the trivial bicomplex 
( [4.75| ) which we used in the context of the sine-Gordon equation. But now we choose a different 
dressing: 

Dz = dz + [5, U]z=-[[1 + 2u x ) z - z] £ + [z y + -(u - u) z] r (4.82) 

where Uz = uz with a field u(x, y). Then we have 5 2 = = <5D + D5 identically, while D 2 = 
is equivalent to 



u xy = -{u - u)(l + 2u x ) . (4.83) 



Since D = d + 5U, following scheme (B) we have R = U2 — TJ\ + T with ST = for a primary 
DBT. The latter condition requires T x = and Tz = Tz. This is satisfied with Tz = az where 
a(y) is a real function. However, since the transformation u% 1— > U\ + a leaves ( |4.83|) invariant, 
we may set a = and obtain 

Rz = {u 2 -u x )z . (4.84) 

The primary DBT conditions now take the form 

(u 2 -u 1 )(l + 2u 2x ) - (u 2 - Ul )(l + 2u lx ) = (4.85) 
(u 2 - ui) y + - (ui + u 2 — U\ — u 2 ) (u 2 - Mi) = . (4.86) 

Adding the first equation to, respectively subtracting it from its complex conjugate, we deduce 

1 1 n U2-UX 1 + (gi + u 2 ) x . s 

«2- u iU = > = --. — —, ; s- ■ (4.87) 

u 2 -ui l + {;ui + u 2 ) x 

Furthermore, one obtains |1 + 2u\ x \ = |1 + 2u 2x \. 

The first permutability condition ( |3.24| ) is identically satisfied and the second permutability 
condition (|3.25|) leads to 

U3 = U * ( U 2-U0)-U! ( Ul -U ) 
U 2 - Mi 

Comparing the operator D with the corresponding operator in the sine-Gordon case, we 
find the transformation 

u*=-(e^-l), <j) y =i{u-u). (4.89) 

Eliminating from the above equations, we obtain ( (4.83j ). If we eliminate u, then we obtain 
the sine-Gordon equation. 

3. Discrete sine-Gordon equation. Let M° be the space of complex functions on an infinite 
plane square lattice and 

(8z) s = (z E - z s )£ + k(z w - z s )t , (dz) s = k(z e - z s )£+ (z w - z s )r (4.90) 
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where k is a real parameter. We use the notation z s = z(x — l,y — l), z E = z(x— 1, y + 1), z w = 
z{x + 1, y — 1), Zn — z(x + 1, y + 1) (see also |T^|). Now d gets dressed with 

(Az) s = [(e^ /2 de-^ /2 )z] 5 = k e^+^V 2 z E £ + e -^-^ )/2 z w r . (4.91) 

The bicomplex condition 5 A = then reads 

e i((f>E-<f>N)/2 _ e i(<l>s-<i>w)/2 — ^2 ^ e «(0 w +0iv)/2 _ & i{<i>E+<l>s)/2\ ^ g2) 

which, multiplied by e l ^ N ~^ E+ ^ w ~ < ^ s ^ 4: , produces the discrete sine-Gordon equation |3^| 

sin[(0iv - <Pe - <Pw + <ps)/4] = k 2 sm[((f) N + (j) E + (j>w + 0s)/4] . (4.93) 
Following scheme (A) of section 3 with az = a z where a is a real constant, we find 



(Rz)s = (92 la 9i z)s =- ,! ' 
Now ( 3.14| ) generates the BT 

' \ + 4>2)e - (01 + 02)5 



*W>l,S+02,s)/2 



sin 



sin 



(02 - 01 V 



— sin 

a 

—a k sin - 



)e + 02 



*>i + 4>2)w + (01 + 02)5 



(4.94) 

(4.95) 
(4.96) 



(see also p2|j ). 

With a w = a 32 = a\ and a 2 o — «3i — «2, (|3.26|) is satisfied and the remaining permutability 
condition (3.25) for the primary DBT reads 



from which we obtain again ( 4.81j) . 



) = «2 (' 



3 i(0O+02)s/2 _ J(<t>l+4>3)s/2 



) 



(4.97) 



^. Infinite Toda lattice. Let M° be the set of real functions Zk(t), fceZ, which are smooth in 
the variable t. On M° we define 



(Sz)k = z k r + (z k+ i -z k )£, (dz) k = (z k - Z k -l) T + z k £ 



(4.98) 



where z = dz/dt. Together with these maps, M = M°®A2 is a trivial bicomplex [I4|. Dressing 
d with A = g _1 dg where g = e Qk , this yields a bicomplex for the nonlinear Toda lattice equation 



(4.99) 



Again, we follow scheme (A) of section 3 to determine a primary DBT. Let g\ = e Pk , g<i = e qk 
and (az) k = a z k _i with a constant a. Then 



- [8(g 2 1 agi)z] k 

a 



9 2 ld (^2^i X ) #1 z 
so that ( |3.14j ) yields 



(p fc _i - g fe ) e"*- 1 " 5 * 2 fc _! r + ( e "+i _ ^-i-ft) (4. 100) 



o'Pk-l-Pk 



-e^-^)z k ^r + (q k -p k )C 



(4.101) 



q k - p k = a ^Pk-ik+i _ gPfc-i-ff*^ 



- g fe = - (e" - e^- 1 "^- 1 ; 



(4.102) 
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We can absorb a in p k by a redefinition p k \— > p k — In \ a\ and choose the sign of t such that the 
above equations become 

q k -p k = _ e Pfc-9fc+i + e Pk-i-g k ) p k i -q k = - e ik- Pk + e <tk-i-Ph-i . (4.103) 

This is a well-known auto-BT of the Toda lattice. From these equations we obtain immediately 

( e <ik-Pky = _ e ik-Qk+i _|_ e Pk-i-Pk r e Ph-i-ihy _ _ e Pk-\-Pk _|_ e ik~±-<ik ^ iq/^ 

Adding these equations and using the Toda equation yields 

( e qk-Pk + tfk-l-toy = _ e 1k-qk+l + e q k -l-q k = ^ (4.105) 

and, after integration, 



q k = e ik-Pk + e p k -i-q k _ lk 



with integration constants 7^. In a similar way, we obtain 



p k = e ik-Pk + e Pk-q k +i _ ^ k 



(4.106) 



(4.107) 



with integration constants 7&. Substituting these expressions into (|4.103| ), we find that 7^ 
7/c = 7 is a constant, and thus 



q k = e qk ~ Vk + e Pk ~ x ~ qk - 7 . 



p k = e qk ~ Pk + e Pk ~ qk+1 - 7 



(4.108) 



which is another form of the auto-BT of the infinite Toda lattice, with a parameter 7 |33 
In terms of gi = e Qi , i = 0, 1, 2, 3, we have 



-qi, k +qj,k-i 



Zk-l 



and with a±o = 0:32 = ct\ and ol 2 q — a u — C12 the permutability conditions amount to 



Q3,k 



-%,k-i + Qi,k + Q2,k + In 



a 2 e 



a\ e' 



Q2,k-l 



a 2 e qi < k — a\ e q2 < k 



(4.109) 



(4.110) 



5. Hirota's difference equation. For functions Zk(u,v) of three discrete variables k,u,v, we 
set (Kz) k (u, v) = z k -i{u,v), (Uz) k (u,v) = z k (u + l,t>), (Vz) k (u, v) = z k (u,v + 1), and define 
bicomplex maps 

6z=(U-1)z£ + (V-1)zt, &z = k x {KU -l)z£ + k 2 {KV -1)zt . (4.111) 
With a dressing similar to that for the Toda lattice we obtain a gauge potential 1-form 

A = e~ q de q = k x (e KU{q) ~ q - 1) KU '£ + k 2 ( e KV{q) - q - 1) KV r (4.112) 
using the notation K(q) = KqK~ x . Now 5 A = with k 1— > k + 1 becomes 



Ki ( 



,q k (u+l,v+l)-q k+1 (u,v+l) _ q k (u+l,v)-q k+1 (u,v) 



) 



«2 (' 



,q k (u+l,v+l)-q k+1 (u+l,v) _ q k (u,v+l)-q k+1 (u,v) 



) (4-113) 
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which, in an equivalent form, is known as Hirota's bilinear difference equation [HJ. 
Following scheme (A), we choose a = aK with a constant a, so that 

R = a e~ q2 Ke qi = a e K{qi) - q2 K . (4.114) 

Now we obtain from ( |3.14j ) the BT 

a ^ e <2l,fcO+M)-92,fc+l(«+l,f) _ e 9l,fe(«,f)-'?2,fc+l(",'")^ 

— Kl ( e Q2,h{ u + 1 ' v )~^,k+l(.u,v) _ e qi, k +l(u+l,v)-q ltk (u,v)^ ^ 115) 
a ^ e <H,k{u,v+l)-q 2)k+1 (u,v+l) _ e qi )k (u,v)-q 2tk+1 {u,v)^ 

— K<2 ^ e q2,k( u ' v+1 )-l2,k+i(u,v) _ e qi, k+1 (u,v+l)-q hk {u,v)^ ^ (4.116) 

As a 'permutability theorem', we obtain the same formula for as in the Toda lattice example. 

6. Principal chiral model. Let M = C°°(1R 2 , C m ) <g> A 2 and 

6z = ztT + z x £, dz = z t T-z x £ (4.117) 

for z G M°. Let G be a group of m x m matrices. Dressing d with the gauge potential 1-form 
A = <7 1 ci <7 , where g G G, this yields a bicomplex for the principal chiral model field equation 

{9~ 1 9x)t + (9~ 1 9t)x = (4-118) 
which is 5 A = (see also ||14j| ). Hence, we follow scheme (A). With 



92 1 &(929i 1 ) 9i = {92 X 92t- 9i l 9it)T - {g 2 1 g 2x - g 1 l 9i x )i (4.119) 
S(92 la 9i) = (g 2 la 9i)tr + (g 2 1 ag 1 ) x ^ , (4.120) 



the primary DBT condition (gig ) and da = requires a to be a constant matrix. Now (|3.14j ) 



results in the following two equations, 

92 X 9it ~ 9f X 9u = (g 2 1 agi)t , 9 2 ~ 1 92 X ~ 9\ X 9ix = -(g 2 la 9i)x • (4.121) 



If a is invertible, the transformation g 2 i— ► a g 2 leaves (|4.118| ) invariant and eliminates a from the 



last equations. This is no longer possible if g is constrained to some subgroup G C GL(m, C). 
For G = U (m) , the last equations read 

02 02* -g\git = {g\agi)t: g\gi x -g\gi x = -{g\ a gi)x ■ (4.122) 

Adding the Hermitian conjugates, the resulting equations can be integrated. Using g\ gi = I, 
they lead to 

glag 1 + g\a^g 2 = C (4.123) 

with a constant real matrix C. Assuming again that a is invertible, it can be written as a 
product u h of a unitary matrix u with a Hermitian matrix h. A redefinition g 2 ^ u g 2 (which 
leaves the field equations and the unitarity constraint invariant) then eliminates u from the 
above equations. Hence we can assume that a is Hermitian. Then we find [g 2 agi,C] = (for 
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all <7i, g 2 ) and thus C = cl with cel. For a = a I with a 6 R, we now recover from ( |4.122j ) 
and ( |4.123| ) a well-known auto-BT for the unitary principal chiral model [33), [35 
With = g^ 1 a,ij gj, the permutability conditions take the following form, 



9s 



(a 3 i gi - a 32 g 2 ) g 1 (g 2 1 



^20 



9i 1 aio] 



-i 



&31 &10 — ^32 &20 



(4.124) 



7. Nonlinear Schrddinger equation. Let M = C°°(R , C ) <8> A 2 with 

^ = z x t + -(cr 3 - 7)2;^, dz = z t r + z x £ (4.125) 

for z G M°. Furthermore, let A = —V r — U £ with 

V = HU. + U>)*, = i (-W fc). (4,26) 

This dressing for d yields a bicomplex for the nonlinear Schrddinger equation 

#f + ^** + 2|V>rV = (4.127) 

(see also [14j]). It is helpful to note that 03 U + U 03 = and U* = —U. The primary DBT 
conditions ( [3.13j ) now imply the following four equations 

R x = V x -V 2 , (4.128) 

l -[a 3 ,R] = U x -U 2 , (4.129) 

R t = V 2 R-RV U (4.130) 

R x = U 2 R-RU 1 (4.131) 

where Ui, Vi are U, V with ip replaced by ipi, i — 1,2. Decomposing R = R + + R~ such that 
o~ 3 R = ±R ± a 3 , we obtain R + = k I + ir 0-3 with functions k and r. Furthermore, ( [4.1 29| ) 
implies 

R~ = i{U 1 -U 2 )a z . (4.132) 
From ( |4.128| ) we obtain k x = and 

tf+ = i(iV> 2 | 2 -|V>i|V 3 . (4.133) 

Hence 

R = R + kI, R = i( . r , (4.134) 

with 

r x = m 2 -m 2 . (4.135) 
From (|4.131|) we obtain k = k, f = r and 

- ^2)* = r (Vi + V2) + « fc - ^2) • (4- 136) 
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With the help of flEI55D , this leads to 



(102 - H% = 2r (|0i| 2 - |0 2 | 2 ) = ~2rr x 

and, after integration, 

r = ±^/a 2 - |0 2 - 0i| 2 
with an integration 'constant' a 2 (t). As a consequence, ( 4.136 ) becomes 

(■01 - 2 ):r = ±(0i + 2 ) a/« 2 - |02 - 0i| 2 + i (-01 - 2 ) 
From (|4.130|) we obtain k t = 0, so that has to be a constant, and 



(4.137) 



(4.138) 



(4.139) 



(0i-0 2 )t = ±i (0i + 2 W" 2 - Vh - 0i| 2 + ~ ^2) (l^i + 2 | 2 + 1^2 - 0i| 2 ) 

-A: (0i - 2 ), • (4.140) 

The last two equations constitute a well-known auto-BT of the nonlinear Schrodinger equation 



27|, |, H, [37|, H-0 

Let us turn to the permutability conditions. First we note that R^ is ant i- Hermit ian, 
traceless and satisfies 



Rij Ru + Rki Rij = [(rij - r M ) 2 + |0j - ipj - fc + 0; | 2 - a 2 - - a 2 kl ] I . 
In particular, we have R 2 , = —a 2 - I- Using trace(-Rjj) = 0, ( |3.24j ) splits into 



R32 
k 3 2 



R31 + -Rio — -R20 

hi + &10 — ho • 



(4.141) 



(4.142) 
(4.143) 



( FT2BD reads 

# 31 .Rio - #32 #20 + feai #10 + fcio #31 - ^32 #20 - ho R32 + hi ho - h2 ho = (4.144) 
and, using ( |4.142| ), 



R3i(Rw — R20 + ho — ho) 

With the help of 



-Rio #20 + (ho~ hi) Rio + (k 3 2~ ho) R 
+ («2 + ^32 ho - hi ho) I ■ 



20 



(Rio - R 20 + ho ~ k 



20 1 



v-1 



-(.Rio - R20 - ho + ho)/p 



where 



p = |0i - 02 1 2 + (r 10 - r 20 ) 2 + {ho ~ ho) 2 



(4.145) 



(4.146) 



(4.147) 



30 The terms proportional to the parameter fc are due to the symmetry transformation ip(t, x) 
e l (- fc2 *+M jp(t, x - 2kt) of the NLS equation. See also @, p. 68. 
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this leads to 



with 



R 31 = - ( a I + b R 10 + c R 20 - 2(k 20 - k w ) R w R 2 o 
P 



(ho ~ k w )[\ipi - ^2 1 2 + (rio - r 20 ) 2 ] + {k 32 - hi) a 



(4.148) 



2 

10 



20 J 



-{ho - k 10 )(k 32 k 20 - k 31 ki + a 
b = -|^i - tp 2 \ 2 - (rio - r 20 ) 2 + « 2 - a 2 + (^20 - k w )(k 

—k 32 k 20 + k 3 i kiQ 
c = al -a 2 w + (k 20 -k w )(k 20 -k 32 ) + k 32 k 20 -k 31 k w . 

Together with R 31 = —a 31 1 this implies 

k 32 = k^ 

and thus k 3 i = k 2Q by use of (KTT43), and moreover 



20 



"31 , 



l 31 — a 20 • 

Now one can derive from ( |4.148| ) the following superposition formula: 



(4.149) 

(4.150) 
(4.151) 



(4.152) 



(4.153) 



^3 



tpo + -{[al - af + [k 20 - k w ) 2 + 2i (k 20 - k w ) r 10 ](^2 - i>o) 



+ [a 2 l0 



« 2 o + (^10 - k 20 ) 2 + 2% (k w - k 20 ) r 2 o](ipi - ipo)} ■ 



(4.154) 



The remaining equations also require a 32 + a 20 = a 31 + af and thus a\ 2 = a 2 . With the 
relations between the parameters of the auto-BT derived above, we recover the permutability 
theorem as formulated in |38 . 



Starting with the trivial NLS solution ip = 0, the BT (fCT3pJ), Q4TH0D determines the 
one-soliton solution 



i[> (x, t) = a j [k x+{a2 " fc2) t+ip] /cosh [a(x — x ) — 2ak t) 



(4.155) 



with real constants cp and xq. Then r = atanhfa (x — x ) — 2a kt]. Let ipj, j = 1,2, be two 
such solutions with parameters atj = cejo,kj = kjo, ipj and Xj (replacing xq). Then the above 
formula for ip 3 determines a two-soliton solution of the NLS equation. 

8. Discrete nonlinear Schrddinger equation. Let M° be the set of C 2 -valued functions Zk(t), 
k G Z, which are smooth in the time variable t. On M° we define 



(5z) k = (z k+1 - z k )r - -(a 3 - I) z k £, (dz) fc = z k r+ (z k - z k ~i)£, (4.156) 

where a dot denotes a time derivative. This determines a trivial bicomplex. Now we dress d to 

(Dz) k = (z k - V k z k ) t + (z k - 2 fc _i - U k z k -i) £ • (4.157) 
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The bicomplex conditions for 5 and D are then equivalent to 



U k+1 -U k --[a 3 ,V k ] = (4.158) 
U k -{V k -V k ^) + U k V k ^-V k U k = 0. (4.159) 

With the decomposition V k = V k + + V k such that <J 3 V k = ±V k a 3 , ( |4.158[ ) implies 

V k =i{U k+1 -U k )a 3 . (4.160) 

In the following we assume that o~ 3 U k = —U k a 3 . Now we decompose ([4.159 ) into ± parts and 
find 

V k + = i U k+1 U k <r 3 , iU k a 3 + (U k+1 + Z7 fc _! - 2U k ) - (U k 2 U k ^ + U k+1 U k 2 ) = . (4.161) 
Assuming furthermore U k ^ = —U k , we can write 



U, 



-tp k 
i>k o 



so that 



V k = i 



-ip k +iip k i[) k+1 -jp k 

1p k+1 - Tpk 1p k+ l1p k 



and from the second of equations ( |4.161| ) we obtain 

i?p k + (V'fc+i + ^fc-i - 2 ^/0 - \i>k\ 2 {i*k+i + V'fc-i) = 



(4.162) 



(4.163) 



(4.164) 



which is the discrete nonlinear Schrddinger equation of Ablowitz and Ladik JIT]. The equations 
( [3.131 ), which determine a primary DBT, with (Rz) k = P k z k _\ lead to 



fc+i 



— Vi,k + Vl,k 

-u 2 , k + u hk 

V2, k P k — P k V\ b-l 



P k — P k -l — U2, k P k -l — P k Uifi-l ■ 

Decomposing P k = P k + P k , ( |4.166| ) implies 

p k = i ( u hk ~ u 2,k) 0-3 

and from Q4.165 ) we obtain 

P k+ i — P k — i (Ui tk+ i Ui tk — U2 tk +i U2,k) o~ 3 . 

Setting 



Vk 
p k 



0-3 



(4.165) 

(4.166) 

(4.167) 
(4.168) 

(4.169) 

(4.170) 

(4.171) 
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we have 

P k = i ( , Pk , (4.172) 

and ( |4.17U| ) becomes 

Pfc+l - Pfc = ^2,fc+l ^2,fc - V>l,Jfc+l V'l.A • (4.173) 

Now we obtain from ( |4.168| ) 

(V>i - ^a)* - Oi - ih)h-\ = Pk-i 4>2,k + Pk ip\,h-\ (4.174) 

and (|4.167|) leads to 

(^1,* - ^2,fc)' = l {lpl,k ~ ?p2,k)(4>l,k4>l,k-l + ^2,fc^2,fe+l) 

+ipk (ip2,k+i ~ ip2,k) + ipk (ipi,k ~ ^i,k-i) ■ (4.175) 

In order to obtain a BT one has to eliminate p k from ( 4.174Q and (|4.175|) with the help of 
( |4.173| ). However, there seems to be no convenient way to achieve that. 

9. A generalized Volterra equation. Let M° be the set of functions z n (t), n G Z, which are 
smooth in the variable t und which have values in C m , m G N. On M° we define 

(5z) n = (z n - z n „ k ) £ + z n r , (dz) n = (z n+1 - z n ) £ + (z n+k+1 - z n ) t (4.176) 

for some fixed k G Z. Then M = M° ® A2 together with 8 and d is a trivial bicomplex. Now 
we introduce a dressing: 

(Dz) n = (g~ 1 dgz) n = (g^ g n+1 z n+1 - z n ) £ + (g~ 1 g n+k+1 z n+k+1 - z n ) r (4.177) 

with an invertible m x m matrix g, depending on t and the discrete variable n. Introducing 
the abbreviation 

V n = g' 1 9n+i , (4-178) 

the operator D can be expressed as follows, 

(Dz) n = (V n z n+1 - z n ) £ + {V n V n+1 ■ ■ ■ V n+k z n+k+1 - z n )r. (4.179) 

The only nontrivial bicomplex condition is <5D + D5 = 0. It results in the generalized Volterra 
equation 

Vn = V n V n+l ■ ■ ■ V n+k - V n . k V n . k+l ■■■V n (4.180) 

which is also known as one of the Bogoyavlenskii lattices ( ||42|| , see also ||43|| ). For k = 1 it 
reduces to (a matrix version of) the Volterra equation 

V n = V n V n+l -V n . 1 V n . (4.181) 
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In order to elaborate primary DBTs, we have to solve the equations fl3.1cp . With the ansatz 
(Rz) n = r n z n+k+ i, we obtain the equations 

r n - r n - k = V 2 , n - V\,n , r n = V 2in ---V 2 ,n+k-Vi, n ---Vi, n+ k (4.182) 

and 

V 2 , n r n +i = r n V hn+k+1 . (4.183) 

Expressing V back in terms of g, we can 'integrate' the last equation and obtain 

r n = 9 2 ,l a 9i,n+k+i (4.184) 

with an arbitrary m x m matrix a(t). Inserted in the two equations Q4.182j ), this leads to the 
following BT for the generalized Volterra equation: 

9 2 ,n 92,n+l ~ 9\,n 9l,n+l = 9 2> n a 9l,n+k+l ~ 9 2 .n-k a 9l,n+l (4. 185) 

fa( g 2,n a 9l,n+k+l) = g2 ) n92,n+k+l— 9i t n9l,n+k+l- (4.186) 

As a permutability relation, we obtain 

93,n = {a 2 gi, n +k+i - ai g 2 , n +k+i) 9o,n+k+i (92,n fl 2 - 9il a iY l (4.187) 

where a\o = a 32 = a\ and 020 = a>z\ = 0,2 and [a%, a 2 ] = 0. 

Let us now restrict our considerations to the scalar case where m — 1. Inserting 

9n = j^- (4.188) 

Jn— fe— 1 



into the BT part (|4.185| ) yields 



fl,n /2,n+l — a fi >n +k+l f2,n-k fl f2,n-k — & fl,n ^2 ,n— 2k— 1 



(4.189) 



fl,n+l f2,n fl,n-k ^2 ,n— k— 1 

which can be 'integrated': 

fl,n /2,n+l ~~ a fl,n+k+l f2,n-k = Pi /l,n+l /2,n (4.190) 

where $ with n = zmod(fc + 1) are 'constants of integration'. The complementary BT part 
( [4.186j ) with a = becomes 

a Dt(fl,n+k+l ' /2,n) — fl,n f2,n+k+l _ a D t {f\, n ' /2,n-fc-l) — fl,n-k-l f2,n ^g-g 
/l,n+fc+l / 2 ,n /l,n /a ,71 — fc— 1 

using Hirota's bilinear operator D t (f ■ h) = fh — fh. This can also be integrated with the result 

a Dt(fl,n+k+l " /2,n) — fl,n /2,71+Jfc+l = 7« /l.n+fc+l /2,7i (4.192) 

where 7« with n = imod(k + 1) are 'constants of integration'. Here we have obtained a BT in 
Hirota's bilinear form. 
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5 Harry Dym equation and equivalence transformations 
of bicomplexes 

Let M = C°°(M 2 , R) <g> A 2 . With 

1 3 

Vz = [(pz x - -<p x z] £ + [3(p 2 z xx + -(tp 2 - 2(p(p xx ) z] r, (5.1) 

Dz = Lp 2 z xx £ + (z t + 4:(p 3 z xxx + 6Lp 2 Lp x z xx )T (5.2) 

we obtain T> 2 = identically, and 

D 2 = ir{-2 V [ Vt + V ^ xxx ]d 2 x } (5.3) 
DD + DD = £r {- [ip t + (f 3 (p xxx ] x - [(fit + ip 3 ip xxx ] d x } (5.4) 

so that the bicomplex conditions are equivalent to the Harry Dym (HD) equation 

(ft + y? 3 Lp xxx = . (5.5) 

A relation between the HD and the KdV equation has been the subject of several pub- 
lications |23[ [PJ. In the following, we show how such a relation emerges in our bicomplex 
framework. This is an instructive example of the application of equivalence transformations to 
bicomplexes. With the gauge transformation 

# = ^-1/22^1/2, D' = v?- 1/2 Dv? 1/2 (5.6) 

(assuming (p to vanish nowhere) we obtain an equivalent bicomplex with 

Vz = ((pd x )z£ + 3(<pd x ) 2 zr (5.7) 
D'z = [(fd x ) 2 + -(2ipip xx - ip 2 )} z£ 

+[d t + A( V d x f + (<pd x ) (2<p<p xx - <p x 2 ) + (p t ]zT. (5.8) 

Next we perform a change of coordinates x — v(s, y), t — s such that v y = <p. As a consequence, 
d y = v y d x = <pd x , d t = d s + y t d y , and = x t = v s + y t v y so that y t = -v s /v y . Writing 
z(t, x) = ((s, y), the bicomplex maps V and D' take the following form in the new coordinates, 

V( = (yZ + 3( yy T, D'( = ((yy-u()Z + (( s + 4( yyy -p( y -q()T (5.9) 

where 

p = Au + ^, q=(lp + Qu) y , u = -l^ + ~(^) 2 = -\s y v (5.10) 



and S y denotes the Schwarzian derivative. The bicomplex conditions for the maps (5.9) reduce 
to 

p = 6u (5-11) 
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where u has to satisfy the KdV equation 



Ug + Uyyy ~ 6 U Uy = . (5-12) 

( pTT| ) together with (pTH|) yields^ 

v s = -(Syv)v y (5.13) 

which must be equivalent to the HD equation. The KdV equation for u is satisfied as a 
consequence of this equation. 

In terms of ip = 1/ \fVy~ the definition of u reads if) yy = uip. Given a HD solution (p(t,x), 
we have to invert y(t,x) = f(l/ip)dx to determine x = v(s,y). Then u = ip yy /ip is a KdV 
solution. Now we can use a KdV-BT to construct a new solution u of the KdV equation. After 
solving tpyy = u ip for ■?/>, we have to invert v = f (l/^ 2 ) dy to find y = y(t, x). Then (p = l/y x is 
again a solution of the HD equation. In particular, if \ satisfies the first of equations ( |4.12| ), i.e., 
Xyy = (u—ct) x with a function a(s), then ip = if) y — (lnx) y ip (which is a Darboux transformation 
|T0| ) satisfies tpy y — uip with u = u — 2(\nx) y y If X a ^ so satisfies the second of equations (|4.12| ) 



(with t, x replaced by s, y), then u is a KdV solution and from if) we obtain a HD solution. 

As an example, let us start with the trivial solution if — 1 of the HD equation. Then we 
obtain y = x + a(t) and thus v = y — a(s) with an 'integration constant' a. Then we have 
ip = 1 and consequently u = which trivially solves the KdV equation. Furthermore, the 
equation Xyy = (u + k 2 )x with a constant k and also the second of equations ( [4.12| ) is solved by 
X = Xo cosh(ky—Ak 3 s). Now we find i/j = — ktanh(ky— 4k 3 s) and-0 = y/k 2 — coth(%— Ak 3 s)/k 3 . 
u = —2k 2 sech 2 (ky — Ak 3 s) is the 1-soliton KdV solution. In order to obtain a HD solution, we 
have to solve x = v{s,y) for y, which results in a function y(t,x) with t = s. This cannot be 
done explicitly, but we find (p = k~~ 2 [l + 1/ sinh 2 (4/c 3 t — ky(t,x))] which indeed solves the HD 
equation. 

Of course, one can try to solve the auto-DBT condition for the bicomplex (M, V, D) associ- 
ated with the HD equation, using ( |1.15| ). This turns out to be rather difficult since already the 



solution for is a non-polynomial differential operator. It appears to be more convenient to 
work with the equivalent bicomplex (M, T>', D'). However, the latter is tied to a less convenient 
form of the HD equation. 



6 Conclusions 

We have introduced the concept of a Darboux-Backlund transformation (DBT) of a bicomplex 
and demonstrated in several examples how BTs for integrable models are easily obtained using 
this simple and universal construction. Once a bicomplex formulation is found for some equa- 
tion, it is straightforward, in general, to apply this method. The bicomplex structure does not 
guarantee a 'decent' BT, however. In some cases, the resulting correspondence between solu- 
tions appears to be practically not of much help (cf the example of the discrete NLS equation 
in section 4.2). 

Higher than primary DBTs have not been sufficiently elaborated in this work, with the 
exception of the Liouville example in section 2. In the KdV case, the secondary DBT turned 

31 This is the simplest case of a Krichever-Novikov equation p5| . 



32 



out to be a composition of two primary DBTs. More precisely, for one choice of sign of a real 
parameter this can only be achieved if one generalizes the primary DBTs to include complex 
transformations. Although the latter do not, in general, generate real solutions from real 
solutions, their composition does. Hence, if we reduce our framework to real solutions and 
real maps, not all of the secondary DBTs are compositions of primary DBTs. Corresponding 
results certainly also hold for higher than secondary DBTs in the KdV case. This shows that, 
in general, we should not expect compositions of primary DBTs to exhaust the hierarchy of 
DBTs. 



Suppose we have three equations EQi 



1,2,3, which are reductions of equations EQ { 



with bicomplex formulations. If is a solution of EQi, then Ui is also a solution of EQ i . Let 
Si and Si denote the solution spaces of EQi and EQ { , respectively. Suppose there are BTs 
BT 21 '■ Si — > <?2 and BT 32 : S2 — > S3 (see Fig. 2) which are determined by primary DBTs. Let 
ui G Si and thus also u\ E Si. Applying BT^BT 21 to it, yields some u 3 G S3 which can be 
projected to some u 3 G S 3 . Not all of such composed and then reduced maps will be trivial 
and not all of them should be expected to be obtainable from primary DBTs of the reduced 
bicomplexes. However, such maps may be recovered as higher than primary DBTs between the 
initial and the final reduced bicomplex. A very simple example is indeed provided by the KdV 
equation mentioned above. In this case, all three equations EQi are given by the real KdV 
equation and EQ i is the complex KdV equation (where the dependent variable has values in 
C). 



Si 



BT 



21 



So 



BT 



32 



7T 2 



713 



< s l 5 2 ► g 3 

Fig. 2: BTs and reductions (the maps 7Tj). 

The method presented in this work is a constructive one, a receipe to determine BTs. We 
do not know, how exhaustive it is and the techniques developed here do not provide us with 
suitable tools to answer this question. Other techniques are available, of course, like those in 
the jet-bundle framework (see 0, chap. 2, and for example). 



In our examples, we have concentrated on equations in two (continuous or discrete) di- 
mensions, with the exception of Hirota's difference equation which depends on three discrete 
variables. Of course, the method also applies to other higher-dimensional equations possessing 
a bicomplex formulation (cf JR], |15j). More examples of this kind will be studied elsewhere. 
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